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Ch ' Abstract 

We study effects of viscous fluid coupled with dark matter in our universe. We consider 
bulk viscosity in the cosmic fluid and we suppose the existence of a coupling between fluid and 
dark matter, in order to reproduce a stable de Sitter universe protected against future-time 
singularities. More general inhomogeneous fluids are studied related to future singularities. 



o . 

O" 1 . I Introduction 

u . 

Recent observational data imply that the current expansion of the universe is accelerating^! . 
This is the so called Dark Energy issue. There exist several descriptions of this fact. Among them, 
the simplest one is the introduction of small positive Cosmological Constant in the framework of 
General Relativity (GR), the so called A-CDM model, but in general also phantom/quintessence 
dark energy (DE) models with effective EoS parameter u approximately equal to —1 are not 
excluded. 

Unfortunately, many of such DE-models bring the future universe evolution to finite-time 
singularity. The classification of the (four) finite-time future singularities has been made in Ref.[2]. 
Some of these types future singularities are softer than other and not all physical quantities (scale 
factor, effective energy density and pressure) necessarly diverge at this finite future time. 

Note that singular solutions correspond to accelerated universe, and often appear as the final 
evolution of unstable de Sitter space. 

The presence of finite-time future singularities may cause serious problems in the black holes 
or stellar astrophysics [3j. Thus, it is of some interest to understand if any natural scenario exists 
in order to cure such singularities. 

In the recent work of Ref.[l], it has been considered the scenario of future singularity removal 
for coupled phantom energy with dark matter and it has been shown the existence of unstable de 
Sitter solution which may solve the coincidence problem but does not cure definitely the singularity. 

The purpose of present work is to generalize the results of Ref.[3], by considering the dark 
matter coupled with viscous fluid, as it is shown in Sec. II. In Sec. Ill we briefly discuss how 
the presence of viscous fluid could produce the singularities and several examples are given (for a 
review, see Ref.[5J). Next, Sec. IV and Sec. V are devoted to reconstruct the standard cosmology 
and the stable de Sitter space avoiding the singularities for this kind of model. In Sec. VI, we 
explore the de Sitter solution in the presence of a generic class of coupled viscous fluids with non 
constant EoS parameter ui and we recover the Cardy-Verlinde formula, whereas in Sec. VII we 
consider several examples of inhomogeneous fluid and the singularity avoidance as a consequence 
of higher-order corrections of modified gravity. Finally, conclusions are given in Sec. VIII. 
We use units of fee = c = ft = 1 and denote the gravitational constant Gn- 
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II Formalism 



Let us recall several fundamental points about GR. The spatially-flat FRW space-time is de- 
scribed by the metric 

ds 2 = -dt 2 + a 2 (t)dx 2 , (III) 

where a(t) is the scale factor of the universe. The Hubble parameter is H — a/a (the point denotes 
the derivative with respect to the cosmic time t). 
The FRW equations of motion (EOM) are: 

3 „ 

H 2 , (II.2) 
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where p and p are, respectively, the total energy density and the total pressure of contents of the 
universe and describe ordinary matter, radiation, dark matter (DM) or more general fluids able 
to explain the accelerated expansion of the universe today. 



From Ea. (|II.2|l and Eq. (|II.3|) . we get the energy conservation law 

p + 3H(p + p) = 0. (II.4) 

It is well-know that ordinary matter and radiation are decoupled and separately satisfy the 
same form of energy conservation law, but it is not necessarily true for other kinds of energy. 

In what follows, we will consider the model where a viscous fluid and dark matter are coupled. 
Their energy conservation laws are given by: 

Pf + 3H(pp + pf) = —Qpf , (II. 5) 

Pdm + iHp DM = Qpf ■ (H.6) 

Here, Q is assumed to be a constant, pdm is the energy density of dark matter (the pressure of 
DM is equal to zero), whereas pp and pf are the energy density and pressure of viscous fluid. 
We assume the fluid equation of state (EoS) in the form[Bl UJ: 

Pf = u{pf)pf - 3H((p F ) , (II. 7) 

3H being the four- velocity of the cosmic fluid. The thermodynamical variable co(pf) is an arbitrary 
function of the density pp- C(pf) is the bulk viscosity and in general it depends on pp- On 
thermodynamical grounds, in order to have the positive sign of the entropy change in an irreversible 
process, ((pf) has to be a positive quantity, so we require C{Pf) > 0. Note that such viscous fluid 
represents special case of more general inhomogeneous fluid introduced in Ref.|13|. 
By taking into account Eq. (|II.7[) . we can rewrite Eq. (|II.5p in the following form: 

Pf + 3H PF (1 + uj(pf)) + Qpf = 9H 2 ({ PF ) , (IL8) 

and, in general, one has C(Pf) as a function of H, ((H). 
The effective EoS parameter uj e f f of the universe is 

where the EOM have been used. If 0J e ff is less than —1/3, d > and the universe expands in 
accelerated way. 

We want to study possible final evolution of the universe due to the presence of dark matter 
coupled with viscous fluid. 

The non viscous case (((H) = 0) has been discussed in Ref.[3], so in this work we suppose 
((H) jLO. 
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Ill Singular future universe 



The presence of viscous fluids in the cosmological scenario could bring the future universe 
evolution to finite-time singularities, in which the Hubble parameter is expressed as: 

H = Tt^W> (IIL1) 

where h and to are positive constants and t < to because it should be for expanding universe. (3 
is a positive constant or a negative non- integer number, so that, when t is close to to, H or some 
derivative of H and therefore the curvature become singular. 

Note that the effect of viscous fluids to Big Rip singularity (f3 — 1 in Eo. piI.l[) ) has been 
discussed also in Ref.[8|. 

Such choice of Hubble parameter corresponds to accelerated universe, because if Eq. pil.l[) is 
a solution of the EOM (|II.2|) - f|II.3|) . it is easy to see that the strong energy condition (p + 3p ^ 0) 
is always violated when /3 > 0, or is violated for small value of t when /? < 0. It means that uJ e ff 
of Ea. pi.9p is less than —1/3 near the singularity. 

The finite-time future singularities can be classified in the following way [5]: 

• Type I (Big Rip): for t — » to, a(t) — » oo, p — » oo and \p\ — > oo. The case in which p and p 
are finite at to is also included. It corresponds to f3 = 1 and f3 > 1. 

• Type II (sudden): for t — > to, a(t) — ¥ ao, p — > po and \p\ — > oo. It corresponds to —1 < ft < 0. 

• Type III: for t — > to, a(t) — > ao, p — > oo and \p\ — > oo. It corresponds to < (3 < 1. 

• Type IV: for t — > to, a(t) — > ao, p — > 0, \p\ — > and higher derivatives of H diverge. The 
case in which p and/or p tend to finite values is also included. It corresponds to j3 < — 1 but 
/3 is not any integer number. 

Many of DE-models may lead to one of such future singularities in the universe evolution. The 
simplest case is the phantom perfect fluid (with equation of state pde = ^>Pde and uj < —1). For 
this kind of fluid, the EOM pl3|) - (|rL3]) admit the solution 

3(1 + LJ) {to - t) 

which corresponds to Big Rip singularity for j3 — 1 and h — —2/3(1 + uj) in Eq. pil.ll) . 

Then, all singularities could appear if we consider the contribute of a non-zero bulk viscosity 
into the fluid equation of state (|II.7|) . 

Let us consider several simple examples[5J, when the thermodinamical parameter of viscous 
fluid uj is a constant, and the fluid is dominant and weakly coupled with dark matter (Q — 0), so 
that its energy conservation law leads: 

p F + 3Hp F (l+uj)~9H 2 (;(p F ). (111.3) 

Suppose to have bulk viscosity proportional to H 

((H)={ H, (III.4) 

where Co is a positive constant. 

The solution of Eq. pil.3[) on the Big Rip singularity (H — h/(to —t)) is: 

PF ~ (2 + 3h(l+uj))(t -t)2 ' { ' 



that solves the EOM when 



H (24nG N (o - 3(1 + uj)) (t - t) ' (IIL6) 
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under the requirement 247rGjvCo — 3(1 + uj) > 0. 

For the same model, also the so called Type I singularity (j3 > 1 in Eo. (|III.H0 could appear. 
In this case an asymptotic solution of Eq. (|III.3j ) is: 

P p ~ 7T~ t — v7 7v28 ' IIL7 

(i + w)(t - typ 

which solves the EOM for u > — 1. 

In the next sections, we consider the model where viscous fluid couples with dark matter and 
we investigate if there exist other realistic possible final scenarios for the universe. 



IV Coupling of viscous fluid with dark matter 

Suppose to have uj constant for the fluid and bulk viscosity in the form 

C - Co# a , (iv.i) 

where a is a real number. In this case, the general solution of Eq. pi.8| ) is 

p -Qt-3u log a(t) q/- p -Qt-3u log a(t) ft /nit') 
e , t "=° e / ji'„Qt +3wloga(t') /.A ■ n'\2 I a \ i I 



PF = POF- 



a(t) 3 a(t) 
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dt'e Qt ' +3ulosa ^a(t')a(t') 2 (^Ty)" . (IV.2) 



where pop is, as usually, a positive constant of integration. 

One possible solution is the de Sitter space, where H = Ho is a constant. One may identify 
the Hubble parameter Hq with the present value of accelerated universe, Hq ~ 10 -33 eV in order 
to reconstruct the standard cosmology. In this case, Ea. (|IV.2[l can be solved as 

(Q + 3H (1 + uj)) 



It follows the solution of Ea. pi.6|) for dark matter 

— n o -3H t n Q „-t(Q+3H (l+u)) , 3ffQ +1 QCo (TV A\ 

Pdm — PoDMe - Pof „ : q jj — e + — — , UV.4j 

where P qdm is a positive constant. It is easy to see that, if Co 7^ 0, the EOM (|II.2|) - ()II.3|) are 
satisfied only if poF = Pqdm = 0. Therefore, we note that, if the de Sitter solution is an attractor 
and is able to describe our universe today, we can require 

Pdm Q 1 /„, r ^ 

— — = TTT ~ o » ( IV - 5 ) 
p F 3H 3 

and the coincidence problem is solved by setting 

Q = H . (IV. 6) 

Now, the ratio of DM and fluid is approximately 1/3, almost independent from initial condi- 
tions. 

By using Eq. pi.2[) - pi.3[) . one has the relation between w and Co necessary to describe the de 
Sitter space: 

w = ~| + |-Ho _1 Co(87rGjv) - (IV.7) 

Here, Eq. (|IV.6[) has been used. With respect to the case of viscosity equal to zero (see Ref.[4]), 
where, in order to solve the coincidence problem u) = —4/3, now the generic condition on uj is 
uj > —4/3, and pf of Eq. pV.3[) is correctly positive. 

For example, a DE-fluid with uj = — 1 admits an exactly de Sitter solution for H = Hq if its 
bulk viscosity is 

J-f a 

C = r , (IV.8) 

and the coupling constant with DM is Q = Hq. 
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V Cosmology 



During the matter dominated era, dark matter has to drive the expansion of the universe. In this 
case, Hq = Q << H. We suppose pp << Pdm, so that Ea. (|II.6[) and Ea. (|II.8[) result decoupled: 

p F + 3H PF (l+uj)~9H a+2 (o, (V.l) 

p D M + 3H PDM * . (V.2) 

The scale factor behaves as 

a(i) = ao t 2/3 , (V.3) 
where ao is a positive constant. In this case, the solution of Ea. (|V.l[) is 

2 xa+2 1 



"* 9 Ms) (1 + 2a ,_ a) ^, (V.4) 



where we have omitted the homogeneus part. Since the solution of Eq. (|V.2|l is 

PODM 

4 

the ratio between energy density of fluid and dark matter when H >> Hq, behaves as 

l-a 



Pdm ^ , (V.5) 



Pf 1 f H 



(V.6) 



Pdm Podm(8ttGn) \ H 

where Eq. pV.7|) has been used. 

It means that, if a < 1, the energy density of fluid increases with respect to the dark matter 
and the universe exits from matter era when the Hubble rate H is close to Hq. 

Note that the condition a < 1 is not strictly necessary. When H becomes close to Hq, the first 
order approximation of Eq. (|IV.2j) on the solution (|V.3P becomes 

2\ Q+2 1/1 1 Hn 



^MsJ (1 + 20,-n) (2 + 2.-n) ^J + ^ ^ 

and so on for the successive orders of Hq: the term of pf proportional to Hq behaves as 1 /t a + 1 ~ n f 
and in any case there exist a value of H close to Hq for which the fluid energy density improves, the 
EOM become inconsistent with respect to the solution (|V.3I) . and the universe exits from matter 
dominated era. 

A possible final scenario is the de Sitter universe described in the previous paragraph (we 
assume that the first term of Eq. pV.4[) is negligible on the de Sitter solution) . 

In order to investigate if the de Sitter solution is an attractor or not, we consider the pertur- 
bation as 

H(t) = H + A(t). (V.8) 



Here, A(t) is assumed to be small. The second EOM (JL3J gives 

2A(t) + 6F A(t) ~ 3H a (a + l)A(t) , (V.9) 
where we have used Eq. (|iV.3[) and Eq. (|IV.7|) . By assuming A(t) — e xt , we find 

A + 3i/ -^ (a + l)~0, (V.10) 

that is 

\~^H (a-l). (V.ll) 

Then, if a < 1, the de Sitter solution is stable and the coupling of viscous fluid and dark matter 
at last generates a stable accelerated universe with a constant rate of DM and fluid. 

If a > 1, the de Sitter solution is not stable and other future scenarios are possible. For 
example, viscous fluid could generate a future singular solution and, if the singularity corresponds 
to the stable solution, the universe could approach to de Sitter space, but will finally evolve to 
such singularity. 
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VI Viscous fluid with non constant omega and the de Sitter 
universe 



Let us consider a more general case, when the thermodinamical parameter uj of viscous fluid is 
not a constant, but it depends on the energy denisty pp, so that u = io(pp). A simple example is: 

w(p F ) = Aopy 1 - 1 , (VI.1) 

where Aq and 7 are constant parameters. The energy conservation law of viscous fluid becomes 

PF + 3HA pl + Q PF = % H a + 2 . (VI.2) 

Here, we suppose the bulk viscosity proportional to H a , ( — (oH a . If we assume 7 >> 1, on the 
de Sitter solution H — Hq, we obtain 

By using Eq. pi.6[) . the energy density of dark matter reads 

PDA! ^ SrPF , (VI.4) 
o-tlo 

and in order to solve the coincidence problem we have to require Q = Hq. 

From the EOM f|II.2[) - 1II.3[) . by assuming the fluid drives the accelerated expansion of the 
universe, it follows 

Jtt +i ( 87rGW x 



A ^3Coi?o tt+i (^7#) . ( VL5 ) 



and by using Eq. (|VI.l[) . one has 

u~ -l + H^CoiSirGN), (VI.6) 

that is an approximation of Eq. pV.7jl . 

In order to investigate if the de Sitter solution is an attractor or not, we consider the pertur- 
bation as in Eq.fVTS]). The second EOM ([TO]) gives 

2A(i) + 6F A(t) ~ H f ) A(t) , (VI.7) 
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where we have used Ea. (|VI.3p and Eq. (|VI.5l) . By assuming A(t) — e At , we find 

\ + 3H -~H (?^\ ~0, (VI.8) 

that is 

^GOt 1 )- 3 )- (vl9) 

Then, if (a + l)/7 < 6, the de Sitter solution is stable. 

It could be of some interests the possibility to recover the Cardy-Verline (CV) formula for this 
kind of fluid on closed universe (with spatial curvature k = 1). A study of the universality of CV 
formula for inhomogeneous fluids has been presented in Ref.[9j (specifically for viscous fluids see 
Ref.[10j). For a perfect fluid with u> constant, the entropy S of the closed universe is written as a 
function of the total energy E and the Casimir energy Ec (for a review, see Ref.|TTj, |12jV 



2na 



^/Ec(2E-E c )) , (VI.IO) 

ap 
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where a is, as usually, the scale factor of the universe and a and /3 are undetermined constants. 
We assume an EoS in the form 

p = A p F -p F - 3Co , (VI.11) 
which is the case as above with bulk viscosity £ = Cq/H. The energy conservation law leads 

a a 

where the prime over p E denotes derivative with respect to the scale factor a. A simple solution 
of this equation is 

''-$)*■ (VL13 » 

The total energy inside the comoving volume V is E = pfV, so that E oc a 3 . For closed universe, 
we can write the energy as the sum of an extensive part Ee and a subextensive part Ec, called 
the Casimir energy, and it takes the form [12 : 

E(S,V) = E E (S,V) + ^E C (S,V). (VI.14) 

The extensive and subextensive parts of total energy, under a rescaling of entropy and volume, 
transform as 

E E (XS, XV) = XE E (S, V) , E C (XS, XV) = X^ 3 E C (S, V) , (VI. 15) 

where A is a constant. 

The FRW universe expands adiabatically, dS = 0, and the products Eca~ 3 and E E a~ 3 should 
be independent on the volume V. Then, by using the rescaling properties of (jVI.151) . we can write 
E E and Ec as functions of the entropy 

E E = ?f, Eo = ^S-*<>. (VI.16) 



By using Ea. (|VI.14l) . we obtain: 



2ira 



-3 x -3 



S= i^=-^/Ec(2E-E c )j • (VI.17) 

This expression is very similar to Ea. (|VI.10[) . The CV formula is recovered for lo = 1, but in 
general only for some special choices or solutions on closed universe, we have an identification 
between the FRW equation of this kind of fluid and the CV formula. 



VII Inhomogeneous dark fluid and singularity avoidance 

In this paragraph we will consider inhomogeneous dark fluid, whose equation of state is given 
by [13] (see also Ref.pj): 

PF = lj{pf)pf + B{ PF , a(t),H, H...) , (VII.l) 

where B is a function of p E , ci(t), H and the derivatives of H. The motivation for this general 
form of time-dependent bulk viscosity comes from the modification of gravity. 

We want to consider several examples, where specific models producing future-time singularity 
are cured by coupled of fluid with dark matter or by adding viscosity term into the EoS. 

Let us start with the following simple EoS: 

pf = -pf + f( PF ) , (vn.2) 
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where J(pf) is a function of the fluid energy density pp. In Ref. j!3| is reconstructed the form of 
/(pf) in order to have the singularity when the fluid energy conservation law is expressed as 

PF + 3Hf{p F ) = 0. (VII.3) 

By using the EOM f|II.2[) - (|II.3|> . we obtain the following solution of the scale factor 

a(t) = a (— ^ , (VII.4) 



to-t 



when the form of J(pf) is 




In Eq. (|VII.4p . n is a positive constant and t < to, so that the Hubble parameter diverges in a 
finite time (t — > to) as in the Big Rip: 

H= nto ( g) 

t(t - 1) V ; 

Therefore, when £ << to, -ff(t) evolves as n/t, which means that the effective EoS parameter 
w e // of Ea. pi.9P is given by w e // ~ — 1 + 2/(3n) > — 1. On the other hand, when t ~ to, 
H(t) ~ n/{to — t) and u> e ff — —1 — 2/(3n) < — 1. The transition between the region w e // > — 1 
and w e // < — 1 occours when t ~ to/2. 

Consider now the coupling of this fluid with dark matter as in Eg. (|II.5[) - (|II.6|) . Since to param- 
eter of DM is equal to zero and violates the condition cj < — 1 for which the Big Rip singularity 
appears, only the fluid could drive the universe to such singularity, so we assume that singular 
solution of Eq. (IVII.6|) already exist and 

3H 2 3n 2 t 2 
PF ~ 8irG N ~ 87rG N t 2 {t - 1) 2 ■ 1 ' 

We note that during the crossing of phantom barrier (w e // = —1), H(t) ~ 2n/t and the 
solution of Eq. pi.6P is: 

pdm — p u . , (vii.8) 

87rGAr(on — l)t 

evaluated in t = to/2. 
We observe that, if 

(6n - 1) t 

the hypothesis that fluid energy density of Eq. (|VII.7p is dominant with respect to the DM energy 
density of Ea. (|VH.8[l when t = to/2 is violated. For large value of coupling constant Q, the uni- 
verse does not cross phantom barrier and is protected against Big Rip future singularity occouring 
when t = to- Note that to is determinated by explicit form of f(p) in Eq. (|VII.5p . 

Let us consider the possibility to obtain the inhomogeneous EoS from the modified gravity. 
The following action is considered: 

S = J d'x^—g + f(R, GO)) , (VII.IO) 

where g is the determinant of the metric tensor and f(R,G) is a function of the Ricci scalar R 
and the Gauss-Bonnet invariant G. By using the variational principle, we can write the FRW 
equations 

H 2 , (VII.11) 



8ttG 



N 
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Here, modified gravity has been included in the energy density p and the pressure p of inhomoge- 
neous fluid so that: 



^ 167T1.7 \ 



^—{24H 3 f G + 6H 2 f' R + 6Hf' R + (f - Rf' R - Gf' G )\ , (VII.13) 
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±— {sH 2 f G + 2f' R + 4Hf' R + l6Hf' G (H + H 2 ) + f R (4H + 6H 2 ) + 

(f-Rf' R -Gf G )}. (VII.14) 

The point denote the derivative with respect to the cosmic time, and we have used the following 
expressions: 



df(R,G) 



dR 



S'g = 



df(R,G) 
dG 



Thus, it follows 
where 

B(H,H.. 



1 



p = up + B(H,H...). 



(l +w )(f-Rf R -Gf G ) 



(VII.15) 
(VII.16) 



167rGjv 

+ {R + !r) \&H 2 {l + w) +4H 

+Hf' R {4 + 6w)+8Hf G \2H + H 2 (2 + 3w) + 2f' R + 8H 2 f G \ , (VII.17) 



and we recover the EoS in the form of Eq. (jVII.ip by writing 

R = 6 f 2H 2 + H) , G = 2AH 2 (H 2 +H 



(VII.18) 



In Ref.|15| are shown several models of f(R, G)-modified gravity producing future-time sin- 
gularities. For example, in the model /(i?, G) = —aG/R, where a is a positive constant, could 
appear the Type I singularity, whereas in the model f{R, G) = ai? 7 , where a and 7 are constants, 
could appear Types II, III or IV singularities. 

It is well know that, in order to cure the singularity problem in f(R, G)-gravity, it is possible 
to use some power function of R, via scenario first suggested in Ref.[16] (for more details, see 
Ref.|17j). or power function of G. As a result, term R m with m > 1 cures Types I, II and III 
singularities. Moreover, if m < 2, the Type IV singularities are cured. The same happens with 
the terms G m : m > 1/2 cures the Type I, II and III singularities, m ^ the Type IV. The Type I 
singularities can be avoided also by using term G m R n , where m and n are positive numbers, and 
so on. Scenario of singularity avoidance for realistic modified gravity is discussed in Ref.[17j, [181 

This kind of terms could be seen as quantum effects in large curvature regime, or could be 
included in the viscosity part B(H, H ...) in the equivalent description of inhomogeneous fluid. 



VIII Conclusion 

In the present paper, we have investigated the possible final scenarios for the universe, in presence 
of viscous fluid coupled with dark matter. 

In principle, viscous fluids could generate any of the known finite-time future singularities, 
which emerge from accelerated universe. If the singularity corresponds to stable solution, the 
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universe could finally evolve to such singularity. A way to prevent this possibility is to consider 
the effects of viscosity on the de Sitter solution descibing the universe today. 

As a result, the presence of bulk viscosity gives rise to a large class of models able to explain the 
current acceleration by producing a stable de Sitter solution protected against singularities. With 
coupling fluid and DM, the coincidence problem could be solved, as first it has been demonstrated 
in Ref.[4], and the viscosity allows to have the stability of De Sitter solution. 

As final remark, we note that inhomogeneous EoS DE/modified gravity may cure singularities 
including corrections in high curvature regime. 
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